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Abstract 

We perform a comprehensive analysis of the decays of charginos and neutralinos in the 
Minimal Super symmetric Standard Model where the neutralino Xi is assumed to be 
the lightest supersymmetric particle. We focus, in particular, on the three-body decays 
of the next -to-lightest neutralino and the lightest chargino into the lightest neutralino 
and fermion-antifermion pairs and include vector boson, Higgs boson and sfermion 
exchange diagrams, where in the latter contribution the full mixing in the third gen- 
eration is included. The radiative corrections to the heavy fermion and SUSY particle 
masses will be also taken into account. We present complete analytical formulae for 
the Dalitz densities and the integrated partial decay widths in the massless fermion 
case, as well as the expressions of the differential decay widths including the masses of 
the final fermions and the polarization of the decaying charginos and neutralinos. We 
then discuss these decay modes, in particular in scenarios where the parameter tan/3 is 
large and in models without universal gaugino masses at the Grand Unification scale, 
where some new decay channels, such as decays into gluinos and qq pairs, open up. 
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1. Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) jl], the lightest neutralinos 
Xi, X2 anc ^ chargino Xi 1 which are mixtures of the higgsinos and gauginos that are the spin 
\ partners of the Higgs and gauge bosons, are expected to be the lightest supersymmetric 
particles. In particular, the neutralino x\ is the lightest SUSY particle (LSP), which because 
of R-parity conservation ||, is stable and invisible. In models where the gaugino masses 
are unified at the Grand Unification scale [||], the masses of these particles are such that: 
m v o ~ m + ~ 2m v o in the case where they are gaugino-like or m^o ~ m + ~ m v o in the 

A2 Xi Al A2 Xi Xi 

case where they are higgsino-like. Thus, the states X2 an d Xi are not much heavier than the 
LSP and might be the first SUSY particles to be discovered. The search for these sparticles 
is a major goal of present and future colliders, and the detailed study of their production 
and decay properties is mandatory in order to reconstruct the SUSY Lagrangian at the low 
energy scale and to derive the structure of the theory at the high scale. 

The decays of charginos and neutralinos have been widely discussed in the literature . 
If the mass splitting between the LSP and the next-to-lightest neutralino X2 or the lightest 
chargino xt is larger than Mz or My/-, the particles will decay into massive gauge bosons 
and the neutralino Xv If n °t) the decays will occur through virtual gauge boson and scalar 
fermion exchanges, leading in the final state to the LSP neutralino and a fermion-antifermion 
pair. Recently, it has been realized [6-9] that for large values of the parameter tan/3, the 
ratio of the vacuum expectation values of the two doublet Higgs fields which are needed to 
break the electroweak symmetry in the MSSM, the Yukawa couplings of third generation 
down-type fermions [b quarks and r leptons] , which are strongly enhanced, lead to dramatic 
consequences for the decays of these particles^. Indeed, the virtual exchanges of, on the 
one side, Higgs particles [because the Higgs boson couplings to b quarks and r leptons are 
proportional to tan /3] and, on the other side, of third generation down-type sfermions [which 
tend to be lighter than the other sfermions in this case] become very important. 

Furthermore, some interest has been recently devoted to models where the gaugino masses 
are not unified at the GUT scale, as it might be the case in a large class of four-dimensional 
string models jTT2| or in the so-called anomaly-mediated SUSY breaking models [I3|. As 



an example, two particular cases have been discussed in Ref. |T3|, where SUSY-breaking 



^^Note that the scenario with tan (3 ~ mt/nrib is favored in models with Yukawa coupling unification at the 
GUT scale [ fLo] . In addition, large tan (3 values, tan/3 > 3-8 depending on the details of the radiative correc- 
tions, are needed to maximize the lightest h boson mass in the MSSM, to cope with the LEP2 experimental 
bound M h > 113.5 GeV [0 in the decoupling regime where the h boson is Standard Model like. 
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occurs via an F-term that is not an SU(5) singlet and in an orbifold string model. In these 
models the gaugino masses at the electroweak scale can be very different from the pattern 
mentioned above. In particular, the X2 an d Xi masses can be closer to the LSP mass in 
some of these models, favoring the occurrence of three-body decays of the light chargino 
and neutralino states [including some some new channels such as X2 ~~ * 9<Z+gluino final 
states], while possibly disfavoring final states with heavy fermions [such as bb final states] 
and therefore affecting dramatically the decay branching ratios. 

In this paper, we perform a detailed investigation of the three-body decay modes of 
charginos and neutralinos in the MSSM, focusing on the scenarios with large values of tan (5 
and with non-unified gaugino masses at the GUT scale. We will provide complete analytical 
formulae for the Dalitz densities of the decays [in terms of the energies of the two final 
state fermions] and for the fully integrated partial decay widths. Furthermore, we will take 
into account the polarization of the decaying particle, which is needed in order to obtain 
the full correlations between the initial state in the production of these particles and the 
final states in their decays. We will also include the dependence on the masses of the final 
state fermions to have a more accurate prediction for final states involving 6-quarks and r 
leptons [especially in scenarios where the mass difference between the decaying particles and 
the LSP is not very large] and to treat properly the case of heavy top quark final states. 
An important ingredient of the analysis will be the inclusion of the effects of the radiative 
corrections to the heavy fermion and chargino/neutralino masses, which will be shown to 
have a large impact. 

This work extends on the recent analyses made in Refs. [6-8] for chargino and neutralino 
decays, and completes our analyses of the higher order decays of SUSY particles [sfermions, 
in particular stops and sbottoms, and gluinos] in the MSSM || |15|| . 

The paper is organized as follows. In the next section, we will summarize the main 
features of the chargino, neutralino, sfermion and Higgs sectors of the MSSM which will 
be needed in our analysis. In section 3, we will display the analytical expressions of the 
(unpolarized) Dalitz densities and the integrated partial three-body decay widths for mass- 
less final state fermions. Section 4 will be devoted to our numerical analysis and a short 
conclusion will be given in section 5. In the Appendix, we present the complete formulae 
for the partial decay widths, including the finite mass of the fermion final states and the 
polarization of the decaying charginos and neutralinos. 
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2. SUSY particles masses and couplings 



To fix our notation, we will summarize in this section the main features of the chargino, 
neutralino, sfermion and Higgs sectors of the MSSM. We will then give, for completeness, 
all the couplings of these SUSY particles [i.e. couplings of the neutralinos and charginos 
to gauge and Higgs bosons and their couplings to fermion-sfermion pairs] as well as the 
couplings of MSSM Higgs and gauge bosons to fermions, which will be needed later when 
evaluating the two-body and three-body partial decay widths. 

2.1 Masses and mixing 

2.1.1 The chargino and neutralino systems 

The general chargino mass matrix, in terms of the wino mass parameter M 2 , the higgsino 
mass parameter \x and tan /?, is given by [PJ 



M 



c 



M 2 y/2M w s p 

y/2M w C/3 /I 



(2.1] 



where we use s@ = sin/5 , cp = cos/5 etc.. It is diagonalized by two real matrices U and V, 

tt* ,i rr i T r I C+ if det.Mc > /n nS 

U*M C V 1 -> U = O- and V = < 1 ., , . . A n (2.2) 

[ a 3 + if detA^c* < K J 

where 03 is the Pauli matrix to make the chargino masses positive and 0± are rotation 
matrices, with angles given by: 

tan 2»_ = 2 t Mw( i' 2C \T 6) • *«» » + = 2 ^<f^ + ^> (2.3) 
This leads to the two chargino masses: 

m 2 xf2 = ~ |M 2 2 + /i 2 + 2M^ T [(Ml - /i 2 ) 2 + AM^(M^c 2 2/3 + Ml + /i 2 + 2M 2f is 2 p)] H(2.4) 

In the limit 3> M 2 , M^, the masses of the two charginos reduce to 

M 2 M 2 
m x ± ~ M 2 - — f (M 2 + /xs 2/3 ) , m x ± ~ + — f e„ (M 2 s 2/3 + /i) (2.5) 

where e M is for the sign of /i. For — > oo, the lightest chargino corresponds to a pure wino 
state with mass m ± ~ M 2 , while the heavier chargino corresponds to a pure higgsino state 

A 1 

with a mass m ± = |w|. 

A2 
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In the case of the neutralinos, the four-dimensional neutralino mass matrix depends on 
the same two mass parameters fi and M 2 , if the GUT relation Mi = |tan 2 #jy, M 2 ~ |M 2 
|T§ is used. In the (-IB, -iW 3 , Hi, # 2 °) basis, it has the form [ C 2 W = I - s 2 w = M^/M 2 Z ] 



M 



N 



Mi -M z s w C{3 M z s w s /3 

M 2 M z c w cp -MzCwS/3 

-M z s w c p M z c w c p -/i 

M Z S W Sp -M Z C W S(3 —\l 



(2.6) 



It can be diagonalized analytically |L7|] by a single real matrix Z . The expressions of the 
masses m Y o are rather involved. In the limit of large lul values, they however simplify to 18 

A 1 



m v o 



*3 



Mi - 
M 2 - 
|/i| + 



M| 
Mf 



(M x + //s 2/3 ) s 



2 



1M| 



(M 2 + /is 2/3 ) c 



2 

w 



2 // 
lMf 
2 /i 



2 



820) // + M 2 s^ + Mic 



2 ^ 



;i + s 2 p) ( 



M 2 s 



Mi c 2 w 



(2.7) 



Again, for — > 00, two neutralinos are pure gaugino states with masses m x o ~ Mi 



xi 

m x o = M 2 , while the two others are pure higgsino states, with masses m x o ~ m x o ~ 
The matrix elements of the diagonalizing matrix, Z^ with i, j = 1, ..4, are given by 



1/4 



Zj 2 

Ail 
Ail 



^2 
Ail 



A^3 



^4 



21-1/2 



{21 



Mi 



tan #w M 2 — ejm v o 



//(Mi - ei m x o)(M 2 - ei m x o) - M 2 s f3 c () {{M 1 - M 2 )c 2 w + M 2 - ei m x o] 

M Z (M 2 - eim x o)s w [iicp + €im x oSp) 
-e i m x o(M 1 - eim x o)(M 2 - e^o) - Mfc 2 /^^ - M 2 )cfy + M 2 - e i m x q 1 



M Z (M 2 - e i m x o)s w [/ic ) g + e^m^s^) 

where is the sign of the ith eigenvalue of the neutralino mass matrix, which in the large 
\/i j limit are: E\ — e 2 = 1 and 64 = — e 3 = e M . Note that we will often use the rotated Zij 
matrix elements: 

Z'a = ZuCw + Z i2 s w , Z' i2 



-ZnSw + Z i2 C\y , ZL — Z a , Z iA — Z,, 



(2.9) 



We will not only discuss the chargino and neutralino spectrum in mSUGRA-type models, 
where the gaugino masses are unified at the GUT scale M GU t, but also when the boundary 
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conditions at this high scale are different. For illustration, we focus on two scenarios discussed 
in Ref. |TJ]]: i) Models in which SUSY breaking occurs via an F-term that is not SU(5) singlet 



but belongs to a representation which appears in the symmetric product of two adjoints: 
(24(g>24) sym =l©24©75©200 [where only model 1 leads to the universal gaugino masses 
discussed previously]. %%) The Oil model which is superstring motivated and where the 
SUSY breaking is moduli-dominated. 

The relation between the gaugino masses at the scale Mguts m i,2,3 ; and at the weak scale 



0(M Z ), M ij2j3 , are approximately given by the relation [19 



M 1 ~ 0A2m 1 , M 2 ~ 0.83m 2 , M 3 ~ 2.6m 3 (2.10) 

leading to the well known hierarchy Mi : M 2 : M 3 = 1:2:6 for a universal gaugino mass at 
the GUT scale, mi = m 2 = m 3 = mi/ 2 , as in mSUGRA type models. The relative gaugino 
masses at Mqut and at the low-energy scale Mz are given in Table 1 ; Ref. fl4| . The pattern 



for the neutralino and chargino masses can be quite different from the universal case 1. In 
particular, for large values of the parameter fi, the LSP is wino-like in the scenario 200 
where M 2 < Mi, implying that Xi and xf are degenerate in mass. In the scenario 75, the 
gauginos x 2 and xf have masses which are very close since |Mi| ~ |M 2 |, while in scenario 
24, the mass splitting between the LSP and the states X21X1 can be very large. In the Oil 
model and if no large loop corrections are present to increase the gluino mass compared to 
the value of M 3 < Mi, M 2 [to avoid the scenario with a gluino LSP], x?> X2 and xf have to 
be higgsino like and can be thus degenerate in mass. 





M 3 


M 2 


Mi 


1 


1(~ 6) 


1(~ 2) 


1(~1) 


24 


2(~ 12) 


-3( 6) 


-1(~-1) 


75 


1(~ 6) 


3(~ 6) 


-5(~ -5) 


200 


1(~ 6) 


2(~ 4) 


10(~ 10) 


Oil 


1(~ 6) 


5(~ 10) 


53/5(~ 53/5) 



Table 1: Relative gaugino masses at Mgut(Mz) in the F$ representations and the Oil model. 

Since X2 and xf c& n be degenerate in mass with the LSP in some of these scenarios, it 
is important to include the radiative corrections to the masses. These corrections are quite 
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involved |20fl . Here we will work in two different approximations, which are valid in the 
(almost) pure gaugino and pure higgsino regions |21], |22|] and which reproduce the complete 
result to better than a few percent. 

For gaugino like neutralinos and charginos, 3> M\,M 2 ,Mz, we will correct only the 
parameters, M 1; M 2 in the chargino and neutralino mass matrices [which means that terms 
of 0(a/A7T x M|//i 2 ) are neglected]; we assume that all fermions are massless and all squarks 
and sleptons are degenerate, with masses nig and mj, respectively; furthermore we work in 
the tree-level decoupling limit for the Higgs sector, where M^ ~ Mz and Mu ~ M H + ~ Ma 
[see section 2.1.3]. For the gluino mass, m g = M 3 + AM3/M3, needed in order to compare 
to the LSP mass, we will include only the dominant QCD corrections. 



In this limit, one then obtains for AM1.2, 3/^1^,3 [^T 

1ST = -£|-{i^.W 2 .°.™^^(M?,o, m ,-)-^ 



S2(3 



(2.11) 



x 



B Q (Ml fi, M A ) - B (Ml fi, M z ) 



B^Ml //, M A ) + B^Ml/i, M z ) 



AM 2 

~m7 



a 



Atts 2 w 



l [ 9B 1 (MlO,m (j ) + 3B 1 (MlO,m I ) - j^s w 



(2.12) 



B (Ml ft, Ma) - B (Ml //, M z ) 



+ B x {Ml n, M A ) + B 1 (Mlfi, M z ) 



-85 (Mf , M 2 , M w ) + AB^Ml M 2 , M w )} 



AM 3 

~mT 



3a s 



{2B (M%, M 3 , 0) - B 1 (Ml M 3 , 0) - 25 1 (M 3 2 , 0, m ? -)} 



(2.13) 



with the finite parts of the Passarino-Veltman two-point functions B\ and B given by 

Log it) -2 



B (q ,m 1 ,m 2 ) 



B l (q 2 ,m 1 ,m 2 ) 



-Log(l 
1 



X 4 



x + Log(l - ./• 

2\ 



-In 



Loefl 



x_Loe:(l 



x 



mj (1 - log^f 



+ (g 2 -m 2 + m 2 )i? (g 2 ,m 1 ,m 2 ) 
with Q 2 denoting the renormalization scale and 



(2.14) 



x± 



2q 2 



q -m 2 + m 1 ±J(q 



m 2 + m 



2 ) 2 -4g 2 (m 2 -ie)) (2.15) 



For higgsino-like x?;X2 an d Xi particles, <C Mi )2 , we will follow the approach of 
Ref. 2^] and only correct the higgsino entries in the neutralino mass matrix and include the 



7 



dominant Yukawa corrections to the light chargino and neutralino masses, due to stop/top 
and sbottom/bottom loopsQ. The masses in the higgsino limit [j22] are then given by [we 
keep the sign of the eigenvalues]: 



777 + ~ 



777^0 ~ 
Al,2 



1 - 



M 2 (u, + 5 C ) 



M 2 



(2.16) 



with 



■in 



X t [Bx((l ,77^, 777^) +Bx(ii ,m t 



>N 



-3a 
~8n~ 



+A 2 (5i(/i 2 , m b , m~ bl ) + 5i(/i 2 , m 6 , m- b2 

A 2 m t s 29t (b (/x 2 , 777,, m h ) - B (u, 2 , m t , m k , 
+X 2 b m b s 2 e b (.BoG" 2 , 7776 , 777^) - B (j2 2 ,m b ,m b2 



(2.17) 



where 0^ are the mixing angles in the stop and sbottom sectors [to be discussed in the next 
subsection] and Xt.b are the reduced Yukawa couplings of the t, b quarks, which in terms of 
the running masses [to be also discussed in the next subsection] are given by: 



A,, 



m b 



y/2MwSwCp 



x t 



y/2M w SwS/3 



[2.18) 



777 v o and 777 v 



m v o are shown 



The xt an d xt masses as well as the mass differences m 
in Fig. 1 as a function of /i for tan/5 = 50, in the five models discussed above. The wino 
mass parameter is fixed to M 2 = 150 GeV and the parameter Mi is obtained from M 2 as 
in Table 1. We see that the mass difference between the lightest chargino and the LSP can 
be very smallQ in models Oil and 200, even after the inclusion of the radiative corrections. 
In model 75, the next-to-lightest neutralino and the lightest chargino can be degenerate in 
mass with the LSP for small values of fi, and the mass difference hardly exceeds 20 GeV [for 
the chosen value of M 2 ] even for large /i values. Note that for values u, > M3, the gluino is 
lighter than the lightest neutralino Xi m model Oil. 

2 We will further approximate the 5c correction by £34 in the paper [[22) . which would be the case for 
almost degenerate squarks; the difference is negligible in general. 

3 The search for charginos and neutralinos, which are almost degenerate in mass with the LSP, can be 
done in e + e~ collisions, either via a search of almost stable particles or by a search of multi-pion final states 
with a large amount of missing energy; see for instance Ref. [Q. At hadron colliders, the direct search of 
such states will be very difficult, if possible at all. 
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Figure 1: The masses of x\ an d Xi an d their mass differences with the LSP Xi as a function 
of n, for tan/3 = 50 and M 2 = 150 GeV with the Mi values given in Table 1. 
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2.1.2 The sfermion system 



The sfermion system is described, in addition to tan j3 and \i, by three parameters for each 
sfermion species: the left- and right-handed soft SUSY breaking scalar masses m^ and 
mj R and the trilinear couplings Af. In the case of the third generation scalar fermions, the 
mixing between left- and right-handed sfermions, which is proportional to the mass of the 
partner fermion, must be included [25]. The sfermion mass matrices read 



M) 



mj + m LL 
m f A f 



m f A f 
ml + m 2 RR 



with ml 



m h = ™% + (4 -efSw) M z C2/3 



RR 



m 



fa 



+ e f s 2 w M 2 z c 2f 3 



(2.19) 



Af = Af — yu(tan/3) 



-21 1 



where // and e/ are the weak isospin and electric charge of the sfermion /, and s 
1 — c 2 ^ = sin 2 9w- They are diagonalized by 2 x 2 rotation matrices of angle Of, which turn 
the current eigenstates, fi and into the mass eigenstates f\ and the mixing angle 
and sfermion masses are then given by 



sin 20 



f 



2m f A f 



m , 



m ■ 



cos 26*4 



m 



LL 



— m 



RR 



2 

m - 



2 

m% 

12 



m\ 
h,2 



m f + - 



in 



LL 



+ m R R T 



m 



LL 



— m 



RR) 



4:mjA 2 f 



(2.20) 



(2.21) 



The mixing is very strong in the stop sector for large values of A t and makes the lightest 
t\ much lighter than the other squarks and possibly even lighter than the top quark itself. 
For large values of tan /3 and \x, the mixing in the sbottom and stau sectors can be also very 
strong, A h T ~ — /itan/3, leading to lighter bi and fi states. 

Since the fermion masses provide one of the main inputs for sfermion mixing, it is impor- 



tant to include the leading radiative corrections to these parameters | 2q| , in particular those 
due to strong interactions. The fermion masses which have to be used in the mass matrices 
eq. ( |2.19|) are the masses rhf{Q 2 ), evaluated in the DR scheme at the scale Q and which, in 
terms of the pole masses mf, are given by pl |: 



m f = i^fiQ 




(2.22) 



In the case of top quarks, it is sufficient to include the one-loop QCD corrections originating 
from standard gluon exchange (first term) and gluino-stop exchange (second term): 
Am t 



m t 



3tt 
3vr L 




+ 5 

B 1 (m § ,m t ~ 1 ) +Bx{m 



(2.23) 



m? 



gi m t 2 > 



s 2 e t — [B (m 
m t 



B (m- g ,m i2 ) 
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where in terms of M = max(mi, 1712), rn = min(mi, m 2 ) and x = m\/m\, the two Passarino- 
Veltman functions [[23| .60,1(^1,^2) = £0,1 (0, rn\,rn%) simply read in this limit 



B (m 1) m 2 ) 
B 1 (m 1 ,m 2 ) 



, (M 2 \ m 2 

log w + l+ ^w 2 



log 



M 2 
Q 2 



log — 
logx 



1 — X 



+ 



x) log X 



(2.24) 



In the case of bottom quarks, the first important correction which has to be included is the 
one due to standard QCD corrections and the running from the scale m b to the high scale 
Q. The DR 6-quark mass [for the NNLO corrections, we assume that the correction in the 



MS and DR schemes are the same, since the latter is not yet available] is given by |27|| : 



rh b (Q 2 ) = rh b (m 2 b ) c[a s (Q 2 )/n} / c[a s (m 2 h )/^] 



with 



m b [m b ) 
c(x) 



5 a s (m 2 b ) 

3 7T 

v 12/23 



12.4 



a 2 (m 2 b ) 



(23x/QY 2/2 - i [l + 1.175x + 1.5x 2 ] for Q 2 < m\ 



c(x) = (7x/2) 4/7 [l + 1.398x + 1.8x 2 ] for Q 2 > m 



(2.25) 



^2.26) 



(2.27) 



After this, one has to include the sbottom-gluino and the stop-chargino corrections which 



are the most important ones [2l|, in particular for large tan/3 and \i values: 
Am b 



"3tt 
a 



m ? , 



Bi(rrig,m bi ) +.B 1 (m s ,mg a ) - s 2 e b — ( B (m s , mjj - B (m~ g ,m b2 

Tfl b \ 

s 26t [BofamfJ - B (ji,m^)] 



Mfc sin 2(3 
M 2 fit&nfl 



a 



A-ks 2 



c 2 9 B (M 2 , m h ) + s 2 e B {M 2 , m h ) + <ji <-> M, 



(2.28) 



V 2 ~ M 2 

For the r lepton mass, the only relevant corrections to be included are those stemming from 
chargino-sneutrino loops, and which simply read 

Am T a M 2 /xtan/3 



.2 2 ^2 [B o (M 2 ,m 0T ) - B (n,mt T ) 



(2.29) 



m T 4irs 2 w /j, 2 - Mf 

The effect of the radiative corrections is shown in Fig. 2 for the case of the bottom quark 
and tau lepton masses for tan/3 = 50 as a function of fi for the various models with and 
without unification of the gaugino masses at Mqut- The wino mass is fixed to M 2 = 150 
GeV and M l5 M 3 ~ m g at the weak scale are given in Table 1. The main correction to the DR 
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bottom quark mass, m{,(Mf) ~ 3 GeV, is due to the SUSY-QCD corrections from gluino- 
sbottom loops in the case of large values of tan (3 and /x. This correction is proportional to 
Arrib ~ — (a s /ir) x tan (3 firrig / m| and can increase or decrease [depending of the sign of /x] 
the 6-quark mass by more than a factor of two. The effect of the radiative corrections is less 
drastic in the case of the r mass since the latter are of the order a few percent. 

Let us now discuss the dependence of the sfermion masses on the gaugino masses as well 
as on the parameters /x and tan /3, in models with a universal mass m for the scalar fermions 
at the scale Mgut ; but without the gaugino mass unification assumption m 12j 3 = mu2- In 
the case of the partners of the light fermions [including 6-quarks], one can neglect to a good 
approximation the effect of the Yukawa couplings in the one-loop Renormalization Group 
evolution of the soft SUSY breaking scalar masses. With the notation of the first generation, 



one then obtains, when including the D-terms, the following expressions [19]: 





= m l 


+ 5.8m 2 3 + 0A7m 2 2 + 4.2 x 10" 3 m 2 + 0.35M| 


cos 2(3 


m k 


= m l 


+ 5.8m 2 + 0A7m 2 2 + 4.2 x l(r 3 m? - 0.42M| 


cos 2(3 


m l R 


= m l 


+ 5.8m 2 3 + 6.6 x KT 2 ™^ + 0.16M| cos 2(3 




m| 


= m l 


+ 5.8m 2 + 1.7 x lQ- 2 m\ - 0.08M| cos 2(3 




m k 


= m l 


+ 0.47m 2 + 3.7 x lO^mj + 0.50M| cos 2(3 






= m l 


+ 0.47m 2 + 3.7 x 10" 2 m 2 - 0.27M| cos 2/3 




m L 


= m l 


+ 0.15m 2 -0.23M§ cos 2(3 





(2.30) 

One has then, in the case of sbottoms and staus, to include the mixing since in this case, 
large enough off-diagonal elements of the mass matrices are obtained for large /x and tan (3 
values [the effect of the trilinear couplings plays only a marginal role]. 

The squark masses are governed by the parameter m.3, while the slepton masses are 
governed by the parameter m 2 , and to a lesser extent m\. Figs. 3a-b show the variation of 
the soft parameters m^ (a) and m ?1 (b) as a function of M 2 for tan (3 = 50 and m = 300 
GeV. As can be seen, depending on the models, the squark and slepton masses can be 
different for different models. In Fig. 3c, the masses mj H ,m e - H and m^^m^ are shown as 
a function of /x for tan/3 = 50; we have used the previous equations and fixed m = 300 
GeV and m\ = m 2 = m 3 = myz = 120 GeV, i.e. as in the mSUGRA-type scenario. While 
for small values of /x, and hence small off-diagonal elements in the b and f mass matrices, 
dR, bi and Br, f\ are almost degenerate in mass, the mass splitting increases with increasing 
/x reaching a substantial amount for /x > m . 
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Figure 2: The b quark and r lepton masses, including the radiative corrections, as a function 
of \i, for tan/3 = 50 and M 2 = 150 GeV in the various models of Table 1. 
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Figure 3: The masses of lightest sbottom (a) and tau slepton (b) as a function of M 2 for 
m = 300 GeV and fi = 750 GeV in the models of Table 1. The b±, dji and f±, e~R masses as 
a function of \x, for tan/3 = 50 and M 2 = 150 GeV in model 1 (c). 
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2.2.3 The Higgs sector 

The MSSM includes two iso-doublets of Higgs fields, which after spontaneous symmetry 
breaking, give rise to a quintet of physical Higgs boson states: h, H, A, H ± f2"8fl. While an 



upper bound of about 130 GeV can be derived on the mass of the light CP-even neutral 



Higgs boson h p9[ , the heavy CP-even and CP-odd neutral Higgs bosons H, A, and the 
charged Higgs bosons may have masses of the order of the electroweak symmetry scale v 
up to about 1 TeV. This extended Higgs system can be described by two parameters at the 
tree level: tan/3 and one mass parameter which is generally identified with the pseudoscalar 
mass Ma- The Higgs mass parameters and the couplings are affected by top and stop loop 
radiative corrections [2jJ, which in the leading approximation are parameterized by 



3Gpmf , m 2 

log^r (2.31] 



V2tt 2 sin 2 (3 m 



where the scale of supersymmetry breaking is characterized by a common squark-mass value 
m. The next -to-leading order QCD corrections can be included by using the running top 
quark mass in the MS scheme. Stop mixing effects can be accounted for by shifting fh 2 in 



eq. (|2.31|) by the amount [A t = A t — fi cot (3] 

rh 2 ^rh 2 + Am 2 : Am 2 = A 2 [l - i 2 /(12m 2 )] (2.32) 

The neutral CP-even and charged Higgs boson masses and the mixing angle a in the neutral 
sector, when expressed in terms of Ma and tan f3, are given in this approximation by 



M 2 A + M 2 z + eT J(M 2 A + M| + e) 2 - 4M 2 Mf c 2 - 4e(M 2 S 2 + Mfc 2 ' 



2/3 _ ^\ 1V1 A*P ' 1v1 Z^Pj 



M 2 H± = M 2 W + M 2 

M\ + Ml TT , 

tan 2a = tan 2(3— - — ^ — —. — with --<a<0 2.33 
H M 2 -M 2 +e/c w 2~ ~ y ' 

In the decoupling limit, Ma 3> Mz, the A, H, H ± bosons become degenerate in mass Ma — 
Mh — Mh± while the lightest h boson reaches its maximal mass value M\ ~ Mf + e; the 
angle a approaches the value a — > (3 — ir/2. The couplings of the h particle to fermions 
and gauge bosons are then SM-like, while the couplings of the H, A, H^ 1 bosons to down 
(up)-type fermions are (inversely) proportional to tan (3. 

In the present analysis, we will use the full renormalization-group improved radiative 
corrections to the Higgs sector given in Ref. [30|. We will often denote the Higgs bosons by 



Hk with k = 1,2, 3, 4, corresponding to H, h, A and H , respectively. 
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2.2 Couplings 

In this subsection, we list the various couplings 0, [16|, 18 1 which will be needed in our 
analysis. All the couplings are normalized to the electric charge e. 

• The couplings of the charginos and neutralinos to the weak gauge bosons W , Z: 



siL,R 



G^ w with 



G R 



-j^[-Z i4 V j2 + y/2Z a Vji] 
l r-[Z l3 U j2 + V2Z i2 U n ) 



(2.34) 



G 



L,R 



Gf-z with 



G L = 1 

ijZ cw$w 

G R - = 1 



~^v i2 v j2 
— \Ui 2 Uj 2 — UnU. 



5ij-s w 



(2.35) 



Gijz with 



i 



(2.36) 



The couplings of charginos and neutralinos to the Higgs bosons: 



G 



L,R 



G 



L,R 
ij4 



with 



C0_ 

s w 

G* =p- 



s w 



ZjaVh + ^ (Z j2 + timOwZji) V i2 



ZjzUn — ^= (Zj 2 + tanOwZji) U{ 



i2 



qL,R 

Xi Xj iJ k 



Xi X-i H 



G ijk With C R 
^ijk 



[ekVixU j2 - d k V i2 Uji} t k 



(2.37) 



r L, R _ r L, R Gi jk = ^{Z j2 -tan6 w Z n ){e k Z l3 + d k Z l4 ) + j 

G X°X°H k ~ G ijk With = _j_ _ t&ndwZji) {ekZi3 + dkZu) £k + • ^ 



where ei j2 = —e 3 = 1 and the coefficients e k and d k read 

ex/di = c a / - s a , e 2 /d 2 = -s a / - c a , e 3 /d 3 = -sp/cp 



(2.39) 



For the couplings between neutralinos, fermions and sfermions, fa — f — X% one nas: 



a 1 
dj 2 

b f 



m f r f { %■ ] _ e f 
V2M w s w \ C9 t J L 

mjrj_{ 



c 9f 



V2M w s w I s e , J R3 X c °s 



(2.40) 
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with r u = Zji/ sin (3 and = Zjs/ cos/3 for up and down-type fermions, and 



ry, $W rj, 

J c w 



(2-41) 



• For the couplings between charginos, fermions and sfermions, fi — f — Xj ■> one nas f° r 
up-type and down-type sfermions: 



ofi \ = }jl S \ + m " y i2 f se u 

b% J " y/2M w s w c fi 1 J 



(2.42) 



+ 



a i2 



1 = m u Vj2 f c 0d 1 



(2.43) 



• Finally, the couplings of the W 7 , Z gauge bosons and the four Higgs bosons H k = H, h, A, H ± 
with k — 1, .., 4 to fermions are: 

/ _2/)-4e / s^ / 2/J ; f 1 

% = ' az = ' ^ = °^ = ( 2M > 

v{= m ^ ,a{ = 0,vl= ^ ,4 = 0,^ = 0, 4= ' m/(ta ^ r2Jj3 (2-45) 
1 2s w M w ' 1 ' 2 2s w M w ' 2 ' 3 ' 3 2s w M w v ; 

j m d t&nf3 + m u cot(3 f m d tan (3 — m u cotj3 
V * ~ 2^[2s w M w ' ° 4 ~~ 2V2s w M w 

with the coefficients r{ 2 as 

ri=s a /sfs , r^ = c a /s p , rf = c a /cf3 , ri = -s a /c p (2.47) 
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3. Three— body decays 



In this section we give the complete analytical expressions of the partial widths of the three- 
body decays of charginos and neutralinos into a neutralino and two fermions, that we will 
denote to be general by u and d [although they can be the same] 

Xi X°j u d (3.1) 

We will not assume that the final neutralino is the LSP x?> but any of the neutralinos x°j 
to cover also the possibility of cascade decays. As shown in Fig. 4, these decays proceed 
through gauge boson exchange [V = W and Z for xt an d X°i decays, respectively], Higgs 
boson exchange [Hk = H + for xt decays and Hk = H, h, A with k = 1, 2, 3 for Xi decays] and 
sfermion exchange in the t- and w-channels [the flavor is fixed by the sfermion-fermion and 
final neutralino vertex]. For gluino decays [ |3"T| , HJ, only the channels with u and t-channel 
squark exchange will be present; the partial widths can be straightforwardly derived from 
those of the neutralino decays, with the appropriate change of the couplings. Note that for 
the treatment of the Majorana nature of the initial state, we use the rules given in Ref. [[£|. 




Figure 4-' The Feynman diagrams contributing to the three-body decays of charginos and 
neutralinos into the LSP and two fermions. 

In this section, we will simply give the complete analytical expressions for the (unpo- 

larized) Dalitz plot density in terms of the energies of two final fermions, and for the fully 

integrated partial widths for vanishing fermion masses^. [In the most complete analysis of 

4 In mSUGRA-type models, this approximation is very good for all light fermion final states, including 
6-quarks and r leptons, since x° an d xt are expected to have masses larger than 0(100 GeV). The ap- 
proximation would be bad for top quark final states; however, if the three-body decays x§ — > x?tt and 
Xi — * Xitb are kinematically allowed, they will not play a major role since the charginos and neutralinos 
will have enough phase space to decay first into the two-body channels \2 ~^ x\Z 1 xlh [ an d possibly XiH 
and XiA] and xt ~ * X^W [ an( i possibly XiH + ], which will be largely dominating. 
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these decays available in the literature up to now, Ref . |J , the fully integrated partial widths 
have not been derived: one integral has been left-out and performed numerically] The for- 
mulae for the general case with non-vanishing values for the masses of the final standard 
fermions [to be able to describe more accurately the cases of chargino decays into rv as well 
as neutralino decays into bb and r + r _ final states and to treat the case of the top quark] 
and where the polarization of the initial gauginos are taken into account, are given in the 
Appendix. 



3.1 The Dalitz densities for the three— body decays 



The Dalitz density of the decay mode eq. (|3.1| ) is given in terms of the reduced energies of 
the two final state fermions 



x 1 = 2E u /m Xi , x 2 = 2E d /m Xi , x 3 = 2E Xj /m Xi = 2 - x x - x 2 
but we will also use the simplifying notation: 

yi = 1 - xi - /i x , y 2 = 1 - x 2 - fJ, x , y 3 = 1 - x 3 + fj, x 



(3.2) 



(3.3) 



with the reduced masses \j? x = M^/m^. [for the final state neutralino we drop the index, i.e. 
// x = m 2 /m 2 }. Neglecting the masses of the final fermions [but not in the couplings] and 



the widths of the exchanged (s)particles, the Dalitz density is given by: 



dr Xi e 4 7?i 



dxidx 2 



64(2vr 



(3.4) 



dr y + dr s + dr d - + dr^ + dr HlH2 + dr Vii + dr yd - + dr sd - + dr$ s + dr $d - 

where N c is the color factor [N c = 3(1) for final state quarks (leptons)] and the dr's corre- 
spond, respectively, to the separate contributions of the square of the gauge boson, u, d and 
Higgs exchanges and the Vu, Vd, ud, <M, $d and H 1 H 2 interferences. 



dr 



The various contributions, in terms of the couplings given in section 2.2, read: 
4 



v 



(v f v - a f v )\G"f + (4 + 4 ) 2 (Gf lV f 



xiyi + 



(4 - a f v ) 2 (Gf tV ) 2 



+(4 + 4) 2 (G^v) 



/^2, 



x 2 y 2 - 4[(4) 2 + (a f v ) 2 }Gf lV Gf iV ^y 3 } (3.5) 



dr,- 



2 

E 



^ (l-s 1 -/i 4 )(l-x 1 -/i J; ) 



(3.6) 
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dr« = E 



^22/2 



^ (l-x 2 - fia k )(l -x 2 - toi) 



(3.7) 



dr* = 2]T 



2/3 



(f{) 2 + (4r 



(V3-fJ>H k ] 



(Gi-jk) 2 + (^fc) 2 P3 + ^\/Jhc{G I j ; jk Gf jk) 



(3.8) 



dr 



Ay 3 v{v{ 



(2/3- PHi)(V3- HH 2 ) 



GijiG ij2 + G {jl G ij2 )x 3 + 2^fJT^{G ijX G ij2 + G^G^) (3.9) 



dr 



yd 



dr 



-4E 



[af fc a? fc Gf iV (4 + 4) + bf k b d k Gf iV (v d - a d v ) 



xiyi 



k=i 



(2/3 -/iy)(l -/ijj 
+ 4) + bf k b%Gf iV (v d - a 



2/3 



4E 

fc=i 



[a^.a" fc G^ y (f 



+ biuWi.Gfj 



V) u ik u jk yj jiV\ V V 



x 2 y 2 



(2/3 - -X 2 - Hu k ) 



TT X [a^ k Gf tV (v^ + al) + bf k b] k Gf lV (vl - a y ) 



2/3 



(3.10) 



(3.11) 



E 



(tfk a u b tk b ji + ^ d ib u jk b d u ){-xiyi - :r 2 2/ 2 + ^32/3) 
(1 -x 2 -A*fiJ(l -^1 - A**) 



+ 



(1-x 2 -A*uJ(1-^i-A*a) J 



(3.12) 



dr w = - E 

k,l 



+ 



dr$ s = E 

k,i 

+ 



' K - 4)<b% (ggfcpEiyi - ^22/2 + 3:32/3) + 2Gf- jky /J^y 3 ) 

(2/3 — A*fc)(l -zi -A*d,) 
>k + { G fjk( x m ~ X 2V2 + x 3 y 3 ) + 2Gf jk ^y 3 



f (vt - aDbltft (Gf^xxyx - x 2 y 2 - x 3 y 3 ) - 2Gf jky /JI^y 3 ) 

I (s/3 — — a?2 — A*«i) 

(at + v^affil (Gf jk (x m - x 2 y 2 - x 3 y 3 ) - 2Gf jk ^y 3 ) \ 



(2/3 — A*fc)(l -X2 — A***,) 
A few remarks need to be made at this stage: 



(3.13) 



(3.14) 
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• In the expressions of the couplings, the indices % and j refer always to the decaying 
chargino or neutralino and the final state neutralino, respectively. 

• For the Higgs boson exchange contributions, in the case of chargino decays, only the 
exchange of the charged Higgs boson is present and in dr$ one has k = 4 only. In the 
case of neutralino decays, the three neutral Higgs bosons will contribute and k in the 
sum J2k °f dr$ runs from k — 1 to 3. In addition, there is an extra term, dT HlH2 , 
due to the interference between the exchange of the two CP-even Higgs bosons h and 
H . Note also that in this case, there is a difference between the contributions of the 
CP-even (and the charged) and CP-odd Higgs bosons which appears in the terms 
£1,2,4 = 1 and e 3 = — 1 in the couplings. 

• For massless final state fermions, there is no interference between the vector boson and 
Higgs boson contributions. In the Appendix, where the fermion mass dependence will 
be included, interference terms between the Higgs bosons and the vector bosons, which 
are proportional to the fermion masses, will be shown explicitly. 

• In the sfermion exchange diagrams, there is a relative minus sign between the ampli- 
tudes of the u and t channels, due to Wick's theorem. This leads to dT V u and dT V( ^ 
contributions which are anti-symmetric in the interchange of X\ and x 2 . In the case of 
dr$„ and dr $ j, the contributions are symmetric in the interchange of x\ and x 2 , due 
to the scalar nature of the Higgs bosons. 

3.2 Integrated three— body partial widths 

Integrating over the energies x\ and x 2 of the two fermions, with boundary conditions, 

l-x 1 -fi x <x 2 <l ^— , 0<X!<l-n x (3.15) 

one obtains the partial decay width, which is given by an expression similar to eq. (3.4): 



a 2 N c 



r v = ra v . 



(3.16) 



r v + Tu + Tj + r$ + t Hi h 2 + + r y j + T-j + r$ s + r 3 

Using the phase space functions and the function Ck defined by: 

\ k = 1 - 2fi x - 2/i fc + (fx k - fi x ) 2 (3.17) 

(3.18) 



Arctan ( ~* +g*T M _ Arctan 
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one has for the various contributions: 



IV = 8 



(4) 2 + (4f (G- V y + (G« v ) 



,A2 



Vx 



6/i 



v 



\ v + fj, v (5 + 5/i x - 7/iy) 



~0- + Ihc - »v)h<Wt - ^(A y + 2/, x )£ y | - 8 [(4) 2 + (4) 2 ] GJ V G* VV ^ 
j 4 0"x -!) + (! + V x ~ 2 Hv)^ogn x + (^Ay - Hv(l +fi x ~ //y)^£y| (3.19) 



I 3 

r/ = E ( a (k4 + b f ikb f u)(°>Wji + Wdl (1 - Hx)(»j k + ^) - 9 (! - A# 
fe,/=i I z 

| (^-1) 2 (^-^ x ) 2 l ^ ^ /; -1 + (^-l) 2 (^/ fc -^x) 2 Log ^/ fc -l 



Log/i x 



(3.20) 



r* = E 2 N) 2 + (4) 2 ][(^) 2 + (G 



r ijk) 



-(1 -/i x )(6/i fe -5-5^ x ) 



4 [-54 



/ife - 3/^ + 7/i 2 + 1 - fi 2 - n x + £ - 5/i fc + 7/J x f4 - 2n x n k 



1-4/ik- 4fi x fi k + 3/4 + /i 2 ) Log/i x + 4 [(v{) 2 + (a{) 2 



j 4 ^ -!) + (! + fhc ~ 2 ^) L °S/ / x + ~ M 1 + A*x - A*k)) A | (3-21) 



H1H2 



2tW{ + [(2 m + 2/^ 2 - 3/i x - 3)(1 - /i x ) 

fl Hl (l + fl x - fl Hl ) fl H2 (l +/i x — AfHa) x r 

i <^Hi*-Hi ^Hnl" 



^H 2 — HHi 
.2 1 ..2 



H 2 ^H 2 



^H 2 ~ 

1+ fl x + fl 2 Hl + H 2 H2 + AiH^Ha - 2(1 + /i x )(^Hx + /!H 2 ))Log/! x 



(/iffi -A*Ha)( 1 +/ i x) 



^H 2 — i^Hx 



^H 2 — VHx 

Log/^ x - 2(1 - ii x ) \ 



\h 2 £h 2 



(3.22) 



Vf 



i^x + 1 - 4 ^/ fc + 2/^y) + 4(1 + _ 2 /i; fc + /iy)A y £ y 
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+ 4 ^x + /4 - 2 ^/ fe - 2 /^ x + 2 Vf k W ~ ^y)Log/i x + (/i x - Hf k ){-1 + Hj k ) 
F{aV, a v _,iXj k ,ii v ) - A f 2 ^/JI^ifi x - 1 - -^Log// x + //y.F(a^, c£, n v ) 

-^(^ - ^ - Wfc + 4)Log-^-) } (3.23) 



^ = A^ k GfA4 + 4) + bt k b%G% v {v d v -a d v )} 

At = < k a] k Gf lV (v u v + a u v ) + blb- k Gf lV (v u v -a^ 

K = -[aia%Gf lV (v d +a d v ) + b d k b%Gf lV (v d -a d v )} 

A u 2 = < k a] k G« v (v»+a») + blb» k Gf iV (vl-a u v ) (3.24) 

\ B( I (1 - /u x )(-l + 2/i fc - fj, x + 2nj { ) - n k \ k C k - fi k (l + n x - fi k )Logfi x 

k,l I V 

-2Log ^' J^ bhc-Vf, ~ IbtVft -2^/,-? r (a? fc ,a-*,^,^fc) j 

-2B{^ U x - 1 - ^Log^ x - _L(^ - ^ - /i x ^ + ^)Log ^ 

V ^/i ^/i ^/i ^x 

+^(a^,a^,^,/x fc )^| (3.25) 



*i = («fc - a u k )a]iblG* k + « + a u k )a^G% k 
B d 2 = (v d k -ai)a d bpf jk + (v d + a d )b d apf jk 

Bl = {vl-al)blalG^ k + {vt + al)alblG? jk (3.26) 



-2 £ {(a>^ + f^x - 1)(2/^ + 2^- - ^ x 

fe,i=i I \ z 



- - !) - ^x L °S^x 



A'wfc A*x A'd, ^"X 



-(/x x - /xs fe /x d - )J r ((fiu k - fi x )/fiu k , fiu k - n x , ji^, jiu k ) \ + (ai k a] k a d ia'ji + b^b^bffi^) 
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U„ 1 /ij 1 

M Lo §^ — ~ ^xX 1 - Vu k )lHu k + Log T7 - i — — - // x )(l - M d ~)//^~ 



A'-Mfc Mx M<2; Mx 

+2(^ x - 1) + (1 + /i x -iiu k )H{Vu k -n x )/Hu k ,Hu k -/i x ,// d -,/i fi J 



In the previous expressions, we have used the variables and functions: 

1 



(3.27) 



(3.28) 



T{a, b, fii, fij) = f(a, Hi) + f(b, Hi) - /(l, /^) + Log^-Log 
F(a, 6, /ij) = f(a, Hi) - f(b, Hi) - f(l, Hi) - Log/^Log 



Hi ~ Hx 
Hi - 1 



/(o,/ii) = Li 2 



Mi - Hx 



Lio 



Mi - 1 



Log(a + fa - l)Log 



Hi- Hx 



(3.29) 
(3.30) 



a + fii — 1J " \a + yUj — 1 y '"' ' ' ' D \/ij-l 
where Li 2 is the Spence function defined by Li 2 (x) = Jq t _1 Log(l — xt)dt. 

3.3 The two— body partial decay widths 

The two-body partial decay widths can be obtained from the expressions given in section 
3.1 by including the total decay widths of the exchanged gauge and Higgs bosons and the 
sfermions. In this case a smooth transition between three- and two-body partial decay 
widths can be obtained. We will list below the integrated form of the two-body partial decay 
widths of charginos and neutralinos into sfermion-fermion pairs [with massive fermions] , and 
into neutralino and gauge or Higgs boson final states; see also Ref. |33[| . 



- ffj) 

+ 



m 



x,. 



(4) 2 + ( & £) 2 )( 1 -^ +/ i /)+ 4 v / M744 aI (m/>/^) ( 3 - 31 ) 



a 



m 



R 

jiV 



(Gf lV ) 2 + (Gf iV ) 2 } (1 + Mx . - fiy) + (1 - Mx . + My )(l - H Xj - Hv)Hv} (3-32) 

1 + Hx 3 ~PH k ) 



a 



Tte - Xj H k ) = ^m Xi xHH Xj ,HH k ){[(Gf jk f + (Gg fc ) 2 
+ ^v^xj GijkGijk} 



(3.33) 



with 



\{x,y) = 1 + x 2 + y 2 -2x -2y -2xy , ^ = 111^/™^ 



(3.34) 
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3.4 Decays into gluino and quark— ant iquark final states 

As discussed in section 2.1, in models without gaugino mass unification at the GUT scale, 
the lightest chargino and the next-to-lightest neutralino could be heavier than the gluino. 
In this case, the three-body decay modes 

Xi -> gud (3.35) 

with Xi = Xi or X2> are kinematically accessible. This decay is mediated by t— and w-channel 
exchange of squarks only; Fig. 5. 




Figure 5: The Feynman diagrams contributing to the three-body decay Xi - > gQQ- 



The Dalitz density and the partial decay width, neglecting the masses of the final state 
quarks, are given by: 



dl\ 



dxidx 2 



e 2 g 2 s m Xi 
8(2tt) 3 



m v . 

4tt Xl 



dr s + dr d - + dr M - 



(3.36) 



with x\ = 2E u /m Xi ,X2 = 2Ed/m Xi . The various amplitudes are as in eqs. ( p^6| , |3~7| , |3.12| ) for 
the Dalitz densities and eqs. ( p.20| , p.27D for the integrated width, with now, // x = m~/m*. 
One has also to replace the final neutralino-/-// couplings, a^bL, by the gluino-quark- 
squark couplings, af, bf, which in the case of mixing read: 



b\ = sin^ , a\ = —b\ = cos9 q 



(3.37) 



The expressions for the three body decays of gluinos into Xi + QQ fi na l states pi], M, are given 



by the previous formulae with the interchange of m Xi and rrig and by dividing the result by a 
factor of 8 to account for the color numbers of the gluino. [Note that this factor is missing in 
the expression of the gluino decay width in Ref . ; since it is a global factor, the branching 
ratios are therefore not affected.!. 
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4. Numerical Analysis 



We will first illustrate our results in an mSUGRA type model, where we assume a universal 
mass m for the scalar fermions and a mass m\/2 for the gauginos at the GUT scale; the soft 
SUSY breaking masses for the Higgs bosons are however disconnected from the one of the 
sfermions so that the pseudo-scalar Higgs boson mass Ma and the higgsino parameter \i are 
free parameters [in contrast to the mSUGRA model where /i is determined, up to its sign, 
from the requirement of electroweak symmetry breaking]. For the squark sector, we will use 
the simple expressions eqs. (|2.30| ) for the soft SUSY breaking left- and right-handed squark 
and slepton masses when performing the RGE evolution to the weak scale at one-loop order 
if the Yukawa couplings in the RGE's are neglected?]. One has then, in the case of the third 
generation sparticles, to include the mixing. Since for sbottoms and stau's, large enough off- 
diagonal elements of the mass matrices are obtained only for large \i and tan (3 values and 
the trilinear couplings play only a marginal role, we will fix the latter to Aj, = A T = —500 
GeV in the entire analysis. We will choose two representative values for tan (3: a "low" value 
(tan /3 — 5) and a large value (tan (3 = 50) and two values for the pseudoscalar A boson 
massQ, Ma = 100 and 500 GeV. 

In a second step, we will relax the gaugino mass unification constraint mi = m 2 = m 3 = 
7711/2 at the GUT scale, and use the weak scale gaugino masses Mi and M 2 given in Table 
1 for the F<s> representations and the Oil model. We will still use the soft-SUSY breaking 
scalar masses given in eq. (|2.30| ). In this case, we will stick in the illustrations to the large 



tan (3 scenario, tan /? = 50, but still show the effect of the Higgs boson contribution by taking 
the two examples Ma = 100 and 500 GeV. 

In most of the cases, the wino mass parameter will be fixed to M 2 = 150 GeV, which 
for large values of /i, leads in an mSUGRA-type model to the masses m x ± ~ m x o ~ 150 
GeV and m x o ~ 75 GeV [there is a very small variation with the value of tan j3] and hence 
to states which are accessible at the high-luminosity phase of the Tevatron and at a future 
e + e~ linear collider with a cm. energy of 500 GeV. 

5 As mentioned previously, for third generation sfermions, neglecting the Yukawa couplings in the RGE is 
a poor approximation since these couplings can be large; this is particularly the case for top squarks which 
however will not be considered in the present analysis, since we will assume that charginos and neutralinos 
are not heavy enough to decay into top quark final states. 

6 In the large tan/3 scenario and in the non-decoupling regime, the experimental bounds on the masses 
of the pseudoscalar Higgs boson A and the lightest Higgs boson h in the MSSM from negative searches at 
LEP2 are M Al M h > 93.5 GeV El. 
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Note that in the entire analysis, we will include the radiative corrections to the 5-quark 
and r-lepton masses, as well as the radiative corrections to the chargino, neutralino and 
gluino masses given in section 2.1. We will also take into account the full dependence on 
the final state fermion masses [using the pole masses mj = 4.6 GeV, m T = 1.78 GeV and 
m c = 1.45 GeV in the phase space, the other fermions are taken to be massless] since in 
some cases [in particular when the decaying chargino or neutralino has a mass which is close 
to the final LSP mass], they play a significant role. 

The branching ratios for the lightest chargino xf and next-to-lightest neutralino x° into 
the LSP and r and 6-quark final states are shown in Figs. 6-8, in model 1 with gaugino 
mass unification at M G ut- The wino mass parameter is fixed to M 2 = 150 GeV and the 
choices tan (3 = 5, 50 and Ma = 100, 500 GeV have been made. 

In Fig. 6a, BR(x ] f — > Xi T+u ) is shown as a function of the lightest fi mass for /i = +500 
GeV. For large values of and with a heavy charged Higgs boson [Ma = 500 GeV leading 
to M H ± = 506 GeV], the branching ratio is small, being at the level of 10%. In this regime, 
the dominant contribution is coming from the virtual W boson exchange and BK(xf) is 
practically the same as BR(W — > //), i.e. ~ 10% for the r + u final state. However, for 
large values of tan (3 and for a light H ± boson [Ma = 100 GeV leading to M H ± ~ 128 GeV], 
the charged Higgs boson contribution [since the H ± ut t couplings are enhanced] becomes 
dominant and the fraction BR(x+ — > Xi T+1J ) can reach the level of 40% even for ~ 500 
GeV. For smaller values of , the virtual stau exchange diagram becomes more and more 
dominant, and BR(x ] f — > Xi T+u ) becomes close to ~ 80% for stau masses of the order of 
150 GeV. If in addition, H ± is relatively light, the branching ratio reaches the level of 100%. 

Figs. 6b and 6c, where BR(%^ — > x°i T+u ) is plotted for a common sfermion mass m = 
300 GeV as a function of /i and tan/5, respectively, show the same trend from a different 
perspective. For small values of tan f3, the mixing in the stau sector and the Yukawa couplings 
of the r lepton are not enhanced and the branching fraction is at the level of 10%. But for 
large tan f3 values, the stau becomes light and the branching ratio becomes close to unity for 
large values of This occurs more quickly, if the charged Higgs boson is light. 

Figs. 7 and 8 show, respectively, the branching ratios BR(%2 ~~ * Xi T+T ~) an d BR(x!] — > 
Xibb), as functions of or m~ bi for /j, = 500 GeV (a), as a function of /i (b) and as a 
function of tan/3 (c) for itlq = 300 GeV. In this case, there is a competition between bb 
and t + t~ final states. In the case of a light A boson and for large tan/3 values, the A 
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and h contributions are much more important in the decay X2 ~^ Xibb than in the channel 
X2 ~^ Xi T+T ~ because of the larger 6-quark mass and the color factor; the Higgs contribution 
makes then BR(x2 - ► X\bb) dominating, except when fi is very light, and the two body decay 
xl — > f\T is close to occur, making BR(x° — ► X\ T+T ~) close to unity. Even for heavy A,H 
bosons, BR(x° — ► Xibb) can reach the level of ~ 50%. However, for large enough values of 
tan/5 and 11, it is the decay channel X2 ~^ Xi T+T ~ which dominates, since for a universal 
scalar mass m , the stau is always lighter than the bi state and its virtual contribution is 
larger, despite of the color factor. Needless to say, the sum of the two branching ratios, 
BR(%2 - *■ Xi T+r ~ + X\bb) 1S m general close to unity. 

In Fig. 9, we illustrate the effect of the radiative corrections to the 5-quark mass [and 
to a lesser extent the tau-lepton mass] by showing the branching ratios BR(%2 ~~ y Xi T+r ~) 
and BR(%2 - ► Xibb) as a function of tan/3 with ii, ttiq, Ma and M 2 fixed to, respectively, the 
values 1 TeV, 300 GeV, 150 and 150 GeV. For li > (< 0), the SUSY radiative corrections 
[in particular, the correction due to sbottom-gluino loops] decrease (increase) substantially 
the value of nib, therefore suppressing (enhancing) the X2 ~^ Xibb rate by a sizeable factor, 
compared to the branching ratio without the correction (solid lines), for large enough tan/3 
values. The fraction BR(%2 — * Xi r+r_ ) increases (decreases) then, accordingly. These 
corrections are therefore very important and must be taken into account. 

In Figs. 10, 11 and 12, we show, respectively, the branching fractions BR(x+ — > Xi t+1> t), 
BR(xa -> XiT + t-) and BR(x° -> x\bb) as functions of //(> 0) in the models 24, 75, 200 and 
Oil without gaugino mass unification as well as in the universal model 1 for comparison. 
The various parameters are fixed to the following values: tan/5 = 50, M 2 = 150 GeV, 
m = 500 and Ma = 100 (a) and 500 GeV (b). Before discussing the various decay channels 
in these models, compared to the universal case, let us make two general comments: 

i) The values of M\^,z at the weak scale are different and modify appreciably the phase 
space for the decays; in particular two-body decay modes and decays into gluinos become 
possible. In addition the radiative corrections to the gaugino masses, although only of the 
order of a few GeV, could allow the opening of channels such as those involving tau leptons. 

ii) Due to the different values of Mi j2 , the evolution of the sfermion masses from A G tjt 
to the weak scale are modified, and the contributions of r sleptons and bottom squarks can 
be enhanced or suppressed compared to the universal case. Also, the radiative corrections 
to the fermion masses are different and can lead to a further enhancement or suppression of 
the Higgs boson and/or sfermion contribution to the decays. 
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Figure 6: The branching ratio BR(xt Xi UT+ ) f or t wo values of 'tan (5 = 5 and 50 and two 
values of Ma = 100 GeV (solid lines) and 500 GeV (dashed lines) as a function of m fl for 
H = 500 GeV (a) as a function of ji assuming m = 300 GeV (b) and as a function of tan (5 
for two values of \i = 100 and 1000 GeV (c); M 2 is fixed to 150 GeV. 
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Figure 7: TTie branching ratio BR(x2 ~^ Xi T+T ) f or t wo values of tan f3 = 5 and 50 and 
two values of Ma = 100 GeV (dashed lines) and 500 GeV (solid lines) as a function of m fl 
for fi = 500 GeV (a) as a function of /i assuming m = 300 GeV (b) and as a function of 
tan/3 for two values of /i — 100 and 1000 GeV (c); M 2 is fixed to 150 GeV. 
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Figure 8: TTie branching ratio BR(x2 ~^ X¥>b) for two values of tan ft = 5 and 50 and two 
values of Ma = 100 GeV (solid lines) and 500 GeV (dashed lines) as a function of m~ bi for 
\i = 500 GeV (a) as a function of /i assuming m = 300 GeV (b) and as a function of tan ft 
for two values of /i — 100 and 1000 GeV (c); M 2 is fixed to 150 GeV. 
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Figure 9: The branching ratios BR(x2 Xi^b (a) and BR(x% — * Xi T+T ) (b) as a function 



of tan 8 for = 1 TeV, M A = 150 GeV, m = 300 GeV and M 2 
without the radiative corrections to the fermion masses. 
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Model 24: For large [i values, \i > 200 GeV, the lightest chargino and neutralinos are 
gaugino like and because M 2 ~ 6 Mi, the mass differences m x + — m x o and m x o — m x o are 
large, making the decays into real gauge bosons, xt ~^ XiW and x\ — > x\Z, kinematically 
possible. The branching ratios for xt an d X2 are then controlled by the W/Z branching 
ratios: BR(W -> r+z/) ~ 10%, BR(Z -> r+r") ~ 3% and BR(Z -> 66) ~ 15%. For smaller 
// values, fj, < 200 GeV, the two neutralinos are mixtures of gauginos and higgsinos and 
three-body decays are possible. The sfermion exchange channels increase the rates for the 
X2 — > Xi T+T ~ an( i Xibb decay channels, with an additional enhancement, in the later channel, 
being to the exchange of the light Higgs bosons for Ma, M h ~ 100 GeV [this contribution is 
milder in the case of t + t~ final states because of the reduced Yukawa coupling]. 

Model 75: For \x ~ (9(200) GeV, m x + — m x o and m x o — m x o are very small even after 
the inclusion of the radiative corrections [Fig. 1] and the decays of X2 an d xt the LSP 
and massive fermions are not kinematically possible [in this case, these particles if they are 
not almost stable, will decay into the LSP and soft pions]. For large values of /i, the mass 
differences between xt, X2 an d the LSP are sizeable [although penalizing the bb final state 
of X2] an< i the different evolution of the sfermion masses as a function of the gaugino masses 
and the different radiative corrections to the bottom and tau lepton masses, explain the 
quantitative differences between the branching ratios in the two models 75 and 1. 

Model 200: Here, the chargino xt an d the LSP are wino-like for large values of //, and 
the mass difference m x + — m x o is too small for the decay xt ~^ Xi T+u to occur. For smaller 
(j, values, this decay can receive large contributions from light Higgs bosons and sizeable ones 
from light sfermions [in particular, b\ is lighter than in model 1]. In the case of the decays of 
the neutralino X21 smce the difference m x o — m x o is always large [exceeding M z for ji > 200 
GeV i.e. when X2 i s bino-like] the branching ratios BR(x2 — ► Xi T+T ~, Xibb) are similar to 
those of model 24 and are controlled by the Z boson decay branching ratios. Note that in 
this scenario, the decay X2 ~ > 9QQ is possible as will be discussed later. 

Model Oil: In this model, the situation is similar to model 200 for the decays of xt- 
Indeed, for \x > 300 GeV, xt is almost degenerate with the LSP and the channel xt ~^ Xi r+U 
is kinematically closed. This is almost the case for the neutralino X2 which has a mass that 
is close to the LSP mass for large /i values, suppressing the bb decay mode. However, the 
new feature in this scenario is that M3 < and for large fi values, the decay modes 

xt,X2iXi ~ * 9QQ open up and become dominant because of the strong interaction part 
[note, however, that the neutralino Xi is n °t the LSP anymore]. 
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Figure 10: The branching ratios BR(xt Xi UT+ ) as a function of /i in models with non- 
universal gaugino masses; we have fixed the parameters to tan/5 = 50, m = 500 GeV, 
M 2 = 150 GeV and M A = 100 (500) GeV for a (b). 
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Figure 11: The branching ratios BR(x2 ~^ Xi r+r ) as a function of n in models with non- 
universal gaugino masses; we have fixed the parameters to tan/5 = 50, m = 500 GeV, 
M 2 = 150 GeV and M A = 100 (500) GeV for a (b). 
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Figure 12: TTie branching ratios BR{x% ~^ Xi^b) as a function of /i in models with non- 
universal gaugino masses; we have fixed the parameters to tan/5 = 50, m = 500 GeV, 
M 2 = 150 GeV and M A = 100 (500) GeV for a (b). 
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Finally, Fig. 13 shows the branching ratio for the decays x\ — * 9QQ i n the model 200 
where m x o < m g < m x o. For small \x values, the lightest neutralinos are higgsino-like and 
they are degenerate in mass. For values of /x around M 2 , the hierarchy m x o < rrig < m x o 
is possible while x?is the LSP, and the decay can occur. However, the neutralino couplings 
to quark-squark pairs are small except in the case of (s)bottoms for large tan/5 values. 
In contrast, the Xi~X2~Z coupling is almost maximal here. BR(x2 gHoq), which is 
approximately the same as BR(%2 - ► gbb), is thus not dominant, but can reach the level of 
25%, despite of the fact that it is a mixed strong-electroweak decay mode. For larger values 
of /i, the neutralinos Xi,2 become gaugino-like and the partial decay widths r(x° —> gqq) 
are more important since the couplings to fermion-sfermion pairs are enhanced; however in 
this case, because M 3 < M 2 , the gluino becomes lighter than the lightest neutralino which 
we assume here to be the LSP. 

In the case of the charginos, the branching ratio for the decays xf ~ * f° r higgsino- 
like charginos is even smaller, since there is no final state with massive fermions [the tb decay 
mode is not kinematically accessible] and the first and second generation (s)particles have 
small couplings for higgsino-like charginos. In the gaugino-like region, the lightest chargino 
becomes lighter than the gluino and the decay does not occur. 
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Figure 13: The branching ratio BR(x\ ~^ 9QQ) as a function of fi in model 200 with non- 
universal gaugino masses. The parameters are fixed to: tan/? = 5 and 50, M 2 = 150 GeV, 
m = M A = 500 GeV. 
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We have developed a fortran code called SDECAY [|35| which calculates the partial decay 
widths and branching ratios of the chargino and neutralino decays. It includes not only the 
three-body decays, X2 ~* Xiff an d Xi ~^ Xiff discussed in this paper, but also all the 
two-body decays of the charginos and neutralinos [including the heavy x° 4 and xt states] 
into gauge bosons, MSSM Higgs bosons and fermion-sfermion pairs. The program contains, 
in addition, the branching ratios for the two-, three- and four-body decay modes of the 
top squarks, as well as the three-body decays of gluinos and all relevant decay modes of 
sfermions other than the top squarks. 

The gaugino mass parameters Mi,M2,Ms, as well as the soft-SUSY breaking scalar 
masses mj L and uif R , can be chosen as free parameters so that decay widths and branching 
ratios can be obtained in non-universal models. However, scenarios with boundary condi- 
tions at high scales are also implemented, since the program has been interfaced with the 



code SUSPECT [36] for the the renormalization group equations for parameter evolution and 
for the proper breaking of the electroweak symmetry. For the parameterization of the MSSM 
Higgs sector, the program has been interfaced with the code HDECAY [37], which in addition 



gives the decay products for the Higgs particles. All radiative corrections discussed in this 
analysis are incorporated into the program. 

We have compared our results with those of Ref. || which have been implemented in 
the program ISA JET |38|| . For massless fermions and if the SUSY radiative correction to the 



fermion masses are not taken into account in the sfermion mass matrices, the agreement was 
very good in models with gaugino mass unification, giving a great confidence that this rather 
involved calculation is correct. [The comparison was slightly involved since the evolution of 
the couplings and the soft SUSY-breaking terms as well as the parameterization of the Higgs 
sectors are given in different approximations in the programs SUSPECT and ISAJET and we 
needed to use the same input parameters at low energy in both programs|].] Our results are 
however different from those which can be obtained with the program SUSYGEN (version 2.2) 
40| used for SUSY particle searches at LEP, since in the latter version, the Higgs boson 



exchange contributions and the effect of third generation sfermion mixing have not been 
implemented^. 



7 We thank Laurent Duflot from ALEPH for his help with this comparison. An independent numerical 
check in the case of the chargino decays into massless final state fermions, has also been performed by F. 
Boudjema and V. Lafage [g9j . 

8 These effects are being included in a new version of the program; we thank S. Katsanevas and N. 
Ghodbane for discussions on this issue. 
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5. Conclusions 



In this paper, we have analyzed the decay modes of charginos and neutralinos in the MSSM 
where the lightest neutralino Xi is the LSP. We focused on the three-body decay modes 
of the lightest charginos xt an d the next-to-lightest neutralinos X2 the LSP and two 
fermion final states, and made a complete calculation of the decay widths and branching 
ratios, taking into account all possible channels: vector boson, Higgs boson and sfermion 
exchange with the mixing in the sfermion sector included. In this context, we have shown 
that the SUSY radiative corrections to the heavy fermion masses, in particular to the b- 
quark mass, and to the chargino and neutralino masses can play an important role. We 
derived full analytical expressions of the Dalitz densities and the integrated partial decay 
widths in the massless fermion case, and provided the complete formulae for the differential 
decay widths, including the finite masses of the final fermions and the polarization of the 
decaying charginos or neutralinos. A fortran code for the numerical evaluation of all the 
branching ratios is made available pi] . 



For large values of tan/?, the bottom and tau Yukawa couplings become large, leading 
to smaller masses of the tau slepton and bottom squark compared to their first and second 
generation partners. At the same time, the Yukawa couplings of tau and bottom quarks 
to the Higgs bosons can become very large. The branching ratios of the decays of the 
lightest chargino into rv final states and of the next-to-lightest neutralino into bb and t + t~ 
pairs can be thus strongly enhanced in this scenario. We have illustrated this possibility in 
mSUGRA-type scenarios where the gaugino masses are unified at the GUT scale, but also 
in scenarios where the boundary conditions for binos and winos are different at this high 
scale, leading to different mass patterns for the charginos and neutralinos, which affect the 
decay branching ratios. In particular, new decay channels, such as the decay of the lightest 
chargino and the next-to-lightest neutralino into gluino and quark-ant iquark final states, 
open up kinematically and can play an important role. 

When SUSY particles will decay via cascades through charginos and the heavier neutrali- 
nos, the events will contain more r leptons and 6-quarks, than first and second generation 
leptons and quarks. This renders the search for SUSY particles and the measurement of the 
SUSY parameters, where the electron and muon channels where used, less straightforward 
as already discussed in Ref. [0]. 6-tagging and the identification of the decays of the tau 
leptons become then a crucial issue in the search and the study of the properties of these 
particles, in particular at hadron colliders such as the Tevatron and LHC. 
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Appendix 



In this Appendix, we will give the lengthy formulae for the three-body partial decay widths 
in the case of finite masses for the fermion final states \p u ^ fid 7^ 0, with fif — m 2 /m 2 ] and 
where the polarization of the decaying chargino or neutralino is taken into account^: 



(A.1) 



where q, p, pi and p<i are the four-momenta of the particles and n Xi is the spin four-vector 
of the decaying "ino" defined by n Xi ■ n Xi = — 1 and n Xi ■ q = 0. 

The partial decay width for both chargino and neutralino three-body decays, following 
the notation given in section 3, is given by: 
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where dTx is decomposed into the spin-independent part [which is half of the unpolarized 
partial decay width] and the part which depends on the spin four vector of the decaying 



"ino": 
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We will use the reduced Mandelstam variables and spin vector: 



u= {q-px) 2 lm x . , 
Spin— independent part: 



t={q~p 2 fjm\ i and n = n x Jm > 



(A.3) 
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[(1 + fi x + fi u + Hd)(u + i) - u 2 - i 2 - 2fi u fi d - fi u - fi d - /i x (2 + fi u + fi d )} 

f f 



+2 \{y s v f - (a f v ) 2 \ ^T d [u + t-fi u - fi d ]j + 2 [(Gf iV ) 2 - (Gf w f 
[u 2 - i 2 - (u - t)(l + fi x + fi u + fid) + {fJ>d ~ - fJ> x )] + 4 G jiV G fivVIbi 
(v f v ) 2 + (a f v ) 2 } [u + t-1- fi x )) - 4 \(v s v ) 2 - (4) 2 1 Vw) } (A.5) 



9 The expressions for three-body decays of charginos and neutralinos, including the polarization of the 
initial states, are available in the literature, see Ref. jf2|, in the case of massless final fermions, no Higgs 
boson exchange and no mixing in the sfcrmion sector. 
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(A.8) 



dr£ = E 



Aa\^/JI^^/fi u fid + 2a^fi x fi u 



k,i=i {~Vd - Hu k + t)(-fi d - fi Ul + f) 
(-fi d - 1 + t) + 2alJfTd(-fiu ~ fi x + *) + <[-^ + *(! + A*x + ^ + A*«) 
-(/x x + /x u )(l + ^)]l (A.9) 



E F7 

fc,Z=l I P% 



- 4aiy/JI^y/tJ, u tJ, d + 2aiy/fi x fi d 



u 



2(4^ffl u ~(- fi d - fi x + u) + a^-u 2 + u(l + fi x + fi d + fiu) 



-(fi x + fl d )(l + flu)] 



(A.10) 
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a l — \ a jk u jl T" u jl u jk 



4 = {4ktfi + 4i hf 3k)(4k4i + HiJ>li) 
4 = (4k4 + b ^i)(a!A + 4b{ k ) 

+ b jk b ji ) K44 + ) ( A - 1 1 ) 



E 77- - i — ^ —rM{jiV)[-{n x + /i,)^ + 1) 

1=1 {1+flx + Hu + Hd-tJ'V-U- t){-jl u - fir + U) { 



-u 2 



+ u(l + fi u + n d + n x )} + b£(jiV)y/jI^(u + i - 1 - fi x ) + b{(jiV)^n u n d 
(u + i- fx u - fi d ) - 4:b{(jiV) v /JI x ~y/fj, u iid + bl(jiV)y/jId~(t - fx x - fx u ) 
+Hti iV )y/tWhc(i - /i d - 1) + 2b f 7 {jiV)^/JI^{u -fi x - (i d ) + 2b{(jiV) % /jr x jT d 

(A.12) 

= E tt- - 4 : r v — {&fO»V) + * - 1 

z=i (1 + A«x + + fid - fiv - u - t)(-fi d - + t) { 

-^x) + b 2(jiV)[-P + t(l + fi u + fi d + fi x ) - (fi x + fi u ){fid + 1)] - 4bi{jiV) 
yffhty/ 'fiufid + b{(jiV)y/(x u (x d (u +i- Hu~ fid) + 2b£(jiV)y/n u fj, x (i - /i d - 1) 
+2b{(jiV)yfjr d {i -thc~ fiu) + b^(jiV)y/fj, x fj, d (u - //„ - 1) + b£(jiV)y/jjZ 

(u-H x -Hd)\ (A. 13) 



b{(ijk) 


= 4A G U< 


+ 4) + 


b W;y 


A 


-4) 


b f 2 (ijk) 


= a{l a jl G ijk( V k 


+ 4) + 




A 


-4) 


b{(ijk) 


= a{i a ji^jk( v l 


-4) + 




A 


+ 4) 


b{(ijk) 


= a{l a j^i J jk( V k 


-4) + 




A 


+ 4) 


b((ijk) 


= 4 b d G U v l + a l) + 


4i b jiG^j^ 


A 


-4) 


b((ijk) 


= 4iHi G U< 


- 4) + 


4i b jiGijk' 


A 


+ 4) 


b f 7 (ijk) 


= a ji b ii G ijk( v l 


-4) + 


4i b jiG^jk' 


A 


+ 4) 


biiijk) 


= 4Mi G U v l 


+ 4) + 


4fijiGf jk > 


A 


-4) 



(A.14) 
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E 7 — J -Tv + ^46*65*] 

+[«>«4 6 tt + - ^ - 1) - 2[a>^ + 

( fl + f _ ^ _ 1) + [a>>^ + blb u ]k b%ai]^jMjT a {i - ^ - 1)1 (A.15) 



dr 17 - 



E 



(f - /i d - 1) + b2(ijk)y/Jj^(-u + fi x + /id) + bj(ijk)y/fj, x fj, d (-u + fx u + l) 

JT d {i - n x - ji u ) + bf(ijk) [u 2 + ui-u(l + fi u + fx d + /i x ) + /i u /i x + ji d ] 



+b c 4 

+26^(ijA;) v //I^ v //i u /id + b*(ijk)y/fi u ti d (u + i- fi u - fi d ) 
+b d 8 {ijk)^(u + t - /i x - 1) I 



(A.16) 



E 



bi{ijk)^fjl^{u - /jl x - fi d ) 



k,l (l+thc + Pu + Pd-V* k -U- t)(-fl d ~ Vui + t) 

+b%(ijk)y/jigi^(-i + 1 + /id) + bl{ijk)y/]T d {-i + /i„ + /i x ) + b^(ijk)y/n x fji d 
(u - 1 - /i u ) + 2b^(ijk)y/Jj^y/fi u fi d + &6(u*0 [«* + ^ - t(l + /i« + /id + /i x ) + A*u 
+/i x /id] + b^(ijk)y/JI^(u + i - n x - 1) + b%(ijk)y/n u (j, d (u + i- /i u - fi d ) \ (A.17) 



Spin— dependent part: 



4) 2 + (4) 



A2 



(1 + /i x + /i n + /id - /iy - W - t) 2 
[pi.n(/i x + /id -m) +p 2 -n(/i x + /i M - t)] +2 (^y) 2 -(af) 2 y/(J, u ii d (pi.n + p 2 .n) 



+2 



( G jiv) 2 + ( G iiy) 2 | ^ y ay[-pi-n(/i x + /^d - w) + P2-n(l*x + /i„ - f)] 



+4Gj iV G^ v v(fa(fy/]I x ~[-p 1 .n(l + ji d - t) + p 2 .n(l + /i u - u) 



(A.18) 
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drf 



= E 



X* (1 + fl x + fi u + fid - - U ~ tf 



vl) 2 + (a{) 2 [l + ^-u-H-2 (v() 



J\2 f„f\* 



(A.19) 



±V f Hl V f H2 



(1 + fx x + fx u + n d - fi Hl - u - + fj, x + fi u + fi d - Hh 2 -u- i) 
G% X G% 2 - G&GgJ [ Pl .n + p 2 .n}(l + fi x - 2^QT d -u-t) (A.20) 



= E 



^jiV^ijk ^jiV^ijk 



k=l (! + fhc + ^ + A*d - V-H k - U ~ t) (1 + H x + fi u + fi d - fly - U - t) 

v k v v + 44} VJ^\Pi- n (i - 1 - /id) + P2-n(-u - 1 
~jT d [p 1 .n{t - fx d + 1) + p 2 .n(-u + n u + 1)] 



+2 



/ / / / 



^jiV^ijk ^jiV^ijk 



44 + 44 



+Pl.ny/fJL d 



f f f f 
K v v ~ a i a v 



(A.21) 



p 2 .n 



^sVJ^K + a ts(t -the- (A.22) 



E 

fej=i 



-pi.n 



{-Vd- V d ,+u)(-n d - n d 



u 



ZatsV^dVx + a ts(u ~Vx~ Vd) \ (A.23) 



where 



4s = ( a jk b ji + 4Mk)(4k4i - Hkrti) 
4s = (4k4i + b jkbii)(-4k4 + b{ k b{i) 



(A.24) 



Kd = E 



i=i (1 + + + - - u - i)(-fi u - fi- di + u) 
+ b 2sU iV )VJ^\Pi- n (i ~ ~ 1) - P2-rc(« 1)] - 2bf s (jiV) v fjI^jI d ~(p 1 .n 

+p 2 .n) + bf s (jiV)y/JI2\pi.n(t - fi d + 1) - p2-ra(« - M« - 1)] 

-2b d %s (jiV)p 2 .n^ fi x fi u - Ab% s (jiV)pi.n^ fi x fi d \ (A.25) 
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E 77- — K \t>WV)y/lhc 

1=1 {I + fJ. x + fiu + fid - fJ-V - U - t)(-Hd - Hu t + t) { 

[pi.ra(f - 1 - - p 2 .n(u - 1 - /i M )] + 2b 2 l s (jiV)p 2 .n(-t + fx x + //„) 
+2b2 s (jiV)^/fx u fx d (p 1 .n + p 2 .n) + Ab^ s (jiV)y/fi u fi x p2.n + 2b^ s (jiV) y /fj, x fj, d p 1 .n 

+K s {jiV)^\p 1 .n{i -fi d -l)- p 2 .n{u -fi u + 1)] 1 (A.26) 



b{ s (ijk) = 4i4 G ?M + a{) - bf^G^ivi - a{) 
b{ s {ijk) = 44 G fM + a{) - b^Gf^l - a{) 

bUijk) = 4Ag«M - 4) - titiPU< + 4) 

b f 4S (ijk) = ada^Gfj^vl - a{) - &£&£G# fc (u£ + a{) 
b{ s {ijk) = a^G'i^ri + a{) - 4Mi G 5k(4 - a{) 

bUm = a^biiG^M - 4) - 4HiGU< + 4) 

b f 7s m = ^b{G% k {v{ - a{) - a^G^ivl + a{) 
b{ s {ijk) = a f ol b{iG? jk (v f k + 4) - 4iHiG- jk {vt - 4) 



(A.27) 



E , - - , * - - J g[^r^W^ 
fc,i=i - A*d, + u)(-fid - Hu k + t) { 

[ Pl .n + p 2 .n] + [a»?^ - %%6gaJ] V^[pi.n(f - ^ + 1) - p 2 .n(£ - ^ - 1)] 

-2[a>^^ - b^b^^T^.n + [a^b^ - a^a^b^^T^.n 

{-t + pL d + 1) - p2.n(-« + ^ u - 1)] + [a^affilb^ - a u lk a%b u jk b d u ] [pi.n(-£ - fi d + 1) 



+p 2 .n(-w - fi u + 1) + [a»^4 - b? k b?Xtfi]y/i^\pi-n(i -f*d-l)-P2 



n[u 



-ix u - 1)] + 2[a£a>J&S - blb^a^^/l^pi.nj (A.28) 

E 7 " =s7 7 1 " 2&is(u'A;) 

fc,i (! + Ibc + + ~ V* k ~ u - t)(-n u - n dl + u) [ 

WxVuP2-n + bjs^ij^y/JIZlpi.nit -fi d -l) + p 2 .n(-u + /J U - 1)] + 2b d zs (ijk) 
y/VxVdPi-n + b^ijty^fJTd^.n^ - /i d + 1) + p 2 .n(-u + fi u + 1)] + b d s (ijk) 
\pi.n(2u + i- 2/i x - - 1) + p 2 .n(u + fi u - 1)] + 2b d s (ijk)y/fi u fi d [pi.n + p 2 .n] 

+b d 8S ( i J k )VK[PiA-i + /id + 1) + p 2 -n(u -fi u - 1)] I (A.29) 
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dr in = Et =s7 ^Ks(y*0 

v//^[pi.n(f - /i d - 1) + P2-n(~u + fJ> u - 1)] - ^b'2 S (ijk) v /fi x fi u p 2 .n + b^ s (ijk) 



y/JId\pi.n(t - fi d + 1) + p 2 -n(-u + /i u + 1)] + 2bl s (ijk)y/fi x fi d p 1 .n + b% s (ijk) 
\pi.n(-i - fi d + 1) + P2-n(-u - 2t + 2/i x + /i u + 1)] + ^(ijfy^fjl^px.n^i + 

/i d + 1) + p 2 -^(w - A*u - 1)] - 2b% s (ijk) v /[i u [id[pi-n + pa.n] | (A.30) 
Phase space 

To obtain the integrated partial widths, one has to express w and t as functions of xi and x 2 

u — 1 - Xi + pL u , t — 1 - x 2 + fid (A.31) 
and integrate over the latter variables, with boundary conditions: 

2^ < Xl < l + [/i u -( v ^+ v ^) 2 ] (A.32) 

s min %2 ^ s max (A. 33) 

1 (Xi - 2)(xi - 1 - fi d + n x - n u ) - \/~K 



2 1 - x 1 + [i u 

1 [X X - 2){X 1 - 1 - \Xd + H x ~ /in) + 

2 1 - X ± + fi u 



with 



A = (/i u - x\) 



(A.34) 



(A.35) 
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